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The impulse response of multimode optical fibers is distorted because 
each mode carries the signal at a different group velocity. Mode coupling 
tends to reduce the width of the impulse response. Rayleigh scattering, 
being the most fundamental scattering process in optical fibers, serves as 
a mode-coupling mechanism. However, it also causes radiation loss. The 
penalty of a seemingly apparent improvement of the impulse response 
through Rayleigh scattering is calculated in this paper. We conclude that, 
because of the high loss penalty, Rayleigh scattering is not a suitable 
technique for pulse-width improvement. 

I. INTRODUCTION 

The term "Rayleigh scattering" describes light scattering from re- 
fractive index inhomogeneities whose linear dimensions are much 
shorter than the wavelength of light. Most of the scattered light 
escapes from the core region of the fiber and enters the cladding or the 
space outside of the fiber. Some of the scattered power goes into other 
guided modes. Rayleigh scattering thus contributes to the losses in 
the fiber and also influences the impulse response through mode 
coupling. 

Since mode coupling tends to improve the impulse response of optical 
fibers, 12 the question may be asked: How beneficial is Rayleigh 
scattering for light transmission in multimode fibers because of its 
mode-coupling ability? To answer this question we investigate the 
loss penalty that is incurred if Rayleigh scattering is assumed as the 
only mode-coupling mechanism. 

For simplicity, our study is limited to a slab waveguide model (see 
Fig. 1) assuming that there is no variation of the refractive index or 
the light field in the y direction. Ignoring coupling between guided 
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Fig. 1 — Schematic of slab waveguide. The scattering centers are distributed 
randomly throughout the core and the outside medium. They are infinitely thin 
threads of slightly different refractive index extending in the y direction. 

modes traveling in opposite directions, we calculate the width of the 
impulse response and the amount of scattering losses. These calcula- 
tions allow us to establish the loss penalty. We find that the loss 
penalty for any significant pulse-width reduction caused by Rayleigh 
scattering is intolerably high. Thus, it is not feasible to improve the 
pulse dispersion of multimode fibers by intentionally implanting 
Rayleigh scatterers into the dielectric material of the fiber. However, 
improved pulse transmission is obtainable by using other carefully 
engineered mode-coupling mechanisms. 2 

II. THE COUPLING COEFFICIENT 

The even guided TE modes of a slab waveguide consisting of a 
perfect dielectric are determined by the y component of its electric 
field. 3 

E v = A cos kx \x\ < d. (1) 

E y = A cos Kde-i^-v \x\ > d. (2) 

The odd guided modes are given by 

E„ = A sin Kd \x\ < d. (3) 



'i/ 



E v = A A sin Kd e" T ( i *'-*> \x\ > d. (4) 



x 



The magnetic field components are obtained by differentiation: 

H,= — '^ (5) 

and 

R y =^- d -^L. (6) 

OJ/Xo ox 

The factor exp [i{ut — /3z)] is omitted from these and all subsequent 
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field equations. The width of the core of the slab is 2d. The parameters 
k and 7 are defined as follows: 

k = (# - /3 2 )* (7) 

and 

7 = (/3 2 - nffc 2 )*, (8) 

with k = a)^[eti^o, n x = core index, and n 2 = cladding index. The 
propagation constant is obtained as a solution of the eigenvalue 
equations : 

tan «d = - for even modes (9) 



K 

and 



tan k(1 = for odd modes. (10) 



7 



The amplitude coefficient is related to the power P carried by the 
mode 

In addition to guided modes, the slab with infinite cladding has radia- 
tion modes. The magnetic fields of the radiation modes follow from B„ 
by means of (5) and (6). The E v component of the even radiation modes 
is* 

E v = B cos ax \x\ < d (12) 

and 

E, = (^jjp) cos[p(|x| -d) +fl \x\>d, (13) 

with \f/ defined by 

a sin ad 
tan^ = — — - -v (14) 

p cos ad v ' 

The amplitude coefficient B is given by 

B / 2 P VoP y 

\ tt/3(p 2 cos 2 ad + (r 2 sin 2 ad) / ' U5J 

The parameters a and p are defined by 

a = ( w ffca _ 0*)i (16) 

and 

P = (n§/c 2 - 2 )*. (17) 

Similarly, for the odd radiation modes we have 

E v = Csin<rx \x\ < d (18) 
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and 

* t -if2ag5y«n|>(i»i-(0 + *] l*l > d - (w) 

Phase is defined by 

Un<t> = t^ (20) 

a COS ad 

and the amplitude factor is given as 

_ / ggw y (21) 

V tt/3(p 2 sin 2 <rd -f <r 2 cos 2 ad) ) 
The coupling coefficient between two modes has the form 5-7 






f° (n 2 - nl)E v Eldx. (22) 

J — 00 



#, and E» are the 2/ components of the electric fields of two modes 
labeled v and m- The index distribution n = n(x, z) describes the wave- 
guide with slight random fluctuations around the average value, and 
n = n (x) is the index distribution that defines the ideal slab wave- 
guide. It is n = ni in the core and n = n% outside. The ensemble 
average of n 2 — no vanishes, 

(n* - nl) = 0. (23) 

The power-coupling coefficients are obtained from the expression 7 - 8 



u-P,)U-z') 



h„ = y f L dz ( L dz'{K.MK*V))* iU, ' r -*- 
L Jo Jo 

= j£4 f dx f dx' f L dz f L dz'({n* - nl){n'> - n'?)) 

1QLP 2 J-co J-*, Jo Jo 

XEX'E'Xe^W^'K (24) 

The prime indicates quantities depending on x' and z'. 

The purpose of this calculation is to study Rayleigh scattering. For 
this reason we may assume that the correlation of the index fluctuations 
reaches only over distances that are much smaller than the wavelength 
2tt/|3,,. The following correlation function is used : 

<(n 2 - nl)(n'> - n« 2 )> = £ 2 <(n 2 - ng) 2 )5(x - x')5(z - z'), (25) 

where D is the correlation length of the index fluctuations. Substitution 
of (25) into (24) leads to 



w 2 eg 



/,,„ ;;iZ) 2 ((n 2 - no 1 ) 2 ) £ \EA*\E»\*dx. (26) 



16P 2 
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The remaining z integration (after integration over the delta function) 
over the distance L resulted in a factor L that canceled from the 
equation. 

To evaluate the remaining integral in (26) we make the following 
assumption. All modes are considered sufficiently far from cutoff so 
that the guided mode fields are very weak at the core boundary, x = d, 
and negligible outside of the core. For a guide supporting very many 
modes, this assumption is justified for most of them. Thus, the integral 
in (26) effectively extends only over the region of the core. The 
integrals are of three different types : 

Ix = / cos 2 K y x cos 2 k m x dx, (27) 

h = / sin 2 k„x sin 2 k^x dx, (28) 

and 

1 3 = I cos 2 k„x sin 2 K^x dx. (29) 

Since almost all modes have rapidly oscillating fields in x direction 
inside of the core, we approximate these integrals by 

Zi-Ii-Ii«g- (30) 

With the help of (11) and (30) we obtain from (26) 

h - - iirrfewi ^ - "° )! >- (31 > 

In the spirit of our approximation, we may assume y„d 2> 1 and 
0, £Zi tiik, where n\ indicates the core index. Thus, the power-coupling 
coefficient can be approximated as follows : 

k 2 
h = h » = 8^d D2{{n2 ~ n2 ° n (32) 

In this far-from-cutoff approximation, the power-coupling coefficient 
is independent of mode number. Rayleigh scattering couples with 
equal strength all of the modes. 

With the same type of approximations, we obtain from (11), (15), 
(26), and (30) the coupling coefficient between a guided mode labeled 
v and an even radiation mode with propagation constant /3 : 

' KPJ 8»i7r 1 1 (p 2 cos 2 ad + a 2 sin 2 ad) K ' 
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Coupling to odd radiation modes leads to the same type of coupling 
coefficient, hl o) (0), except that cos ad and sin ad are now interchanged. 
The power (scattering) loss coefficient for mode v is 



f njfc 

a. = 2/ LhPW+hPWy* (34) 

Jo 

This expression can be justified as follows. The power-coupling co- 
efficient indicates the amount of power flowing per unit length from 
the guided mode to each individual radiation mode. The sum of the 
contributions to all radiation modes gives the total loss. Since radia- 
tion modes form a continuum, the sum becomes an integral. The 
factor 2 in front of the integral indicates the doubling of the loss 
caused by power flowing not only into forward but also into backward 
traveling radiation modes. The integral over p can be converted to 
integration over /3 as follows : 

«. = 2 /" "'* W 08) + W 08) ] - dp. (35) 

Jo p 

The integration includes only propagating radiation modes. The con- 
tribution of even and odd modes is very nearly the same, so that we 
use only the coupling coefficient (33) and double the factor in front 
of the integral : 

TO 2 ((n 2 - ng) 2 ) /•"'* P d0 m) 

a " 2vn l Jo p 2 cos 2 ad + a 2 sin 2 ad K ' 

To the approximation used in this analysis, the power-radiation-loss 
coefficient of the guided modes is independent of mode number. 

An exact solution of the integral in (36) is hard to obtain. If we 
consider the fact that for large values of d the sine and cosine functions 
pass through many periods throughout the range of integration, we 
can replace the integrand by its average value over a few periods 
of the periodic functions. This average is 

/ P \ m 1. (37) 

\ p 2 COS 2 ad + a 2 Sin 2 ad /average a 

It now remains to solve the integral : 

f n » k dP f nik dP . 7l 2 /qOX 

/ — = I , = = arcsin — , (38) 

Jo <r Jo Vn?fc 2 - /3 2 n i 

In most cases of practical interest the ratio n 2 /n l is very close to unity 
so that we can approximate the integral by *-/2. We thus obtain the 
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following equation for the radiation power loss coefficient 



a = a„ = 



= 2nikhd. (39) 



4ni 
The last part of the equation follows from (32). 

III. CALCULATION OF IMPULSE RESPONSE 

Pulse propagation in optical fibers can be described by the following 
equation for the average power 2 - 7 

dP„ 1 dP» N 

This system of coupled power equations holds only for modes traveling 
in the same direction. Rayleigh scattering scatters power in forward 
as well as backward directions; however, we must ignore the back- 
ward scattered power flowing into guided modes. Physically, it appears 
that this approximation should pose no difficulty, since only those 
modes that travel in near synchronism have a chance to interact 
thoroughly. Backward scattered power travels away from the pulse 
that created it; thus, it cannot alter the shape of the impulse response 
except, perhaps, by repeated reflections. Backward scattered power 
contributes mainly to the scattering losses. We have taken backward 
scattering into radiation modes into account, but the additional loss 
caused by backward scattering into guided modes contributes far less 
loss and is ignored in our treatment. Thus, we recognize that the ap- 
proximation may lead to a slight underestimation of the total scatter- 
ing loss. 

To solve (40) we use the trial solution 

P v = A ,*-»*•*- i*-<»i-/rfl . (41) 

Substitution into (40) leads to 



A > = h {k A )/[«-° + Nh + i "{v,- n i)\ 



(42) 



We used the fact that the loss coefficients and the coupling coefficients 
are independent of the mode number. The quantity N is the total 
number of modes. 

We obtain the group velocity of the modes from an approximation of 
the propagation constant. Using 9 

k4 « v \ , (43) 
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we obtain from (7) 

'*[*-('h)T" (44) 

The inverse group velocity of mode v is 

1 -= ^ = **§. = ^. (45) 

y„ do) c dk cf3 

Using niA; 2> cir/2d we obtain approximately 



with 



I «?*(1 +(?,»), (46) 

u„ c 



G = oT^TwvT < 47 > 

8{nikd) 2 



The solution of the equation system (42) is accomplished with ease, 
since we realize that the sum term in the numerator is independent of 
the mode label. Thus, the coefficients A, must be of the form 

A r = - (48) 

a - a + Nh + io> - Gv 2 
c 

Substitution of (48) into (42) leads to an eigenvalue equation for the 
determination of <r : 

h E — = 1- (49) 

H " l a - a + Nh + w - Gn 2 
c 

The sum can be approximated by the integral 

r N dx = 1 

Jo tt _ r + jyfc + fa^« |"iw^(?(o-<r + ^)J i 

VI 

(50) 




X arctan v , v . 

\a. — a + i\n 

Thus, we obtain from (49) and (50) the eigenvalue equation 



iw^GN 2 Y 



Fortunately, we need only the lowest-order eigenvalue since it has the 
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significance of the steady-state loss of the system of coupled modes 
and also determines the shape of the impulse-response function. 2 - 7 
The solution of (51) is accomplished by using the fact that the lowest- 
order eigenvalue must be close to the loss coefficient a. Thus, we set 

tr = a + n. (52) 

Next, we expand the tangent function in series and solve for t\. In 
this way we obtain the approximate solution 

For our purposes the coefficient p of o> 2 is of most importance. We 
obtain the general pulse shape by substituting (53) into (41) and 
integrating over w from — « to «> . Neglecting the a> dependence of 
A„ we find a Gaussian-shaped pulse whose width is 27 

The width of the signal in the absence of mode coupling is 

AT = (— - -)l = ^ GN*L. (55) 

\V N ViJ c K ' 

The relative improvement of the width of the steady-state pulse in 
the presence of mode coupling is expressed by the factor 2 ' 7 

7? - M - 8 KM 

Using (39) and (56), we define the loss penalty by the expression 2 - 7 

K*aL = 2.8 ^P- (57) 

The number of modes is obtained from (44) with the help of the cutoff 
condition /3 = n 2 k for v = N ; thus, 

A^ = — V^T^I. (58) 

IT 



The expression for the loss penalty thus assumes the form 

*•* ~ 3.1 



4 4 31 
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IV. DISCUSSION 

We can now answer the question that was asked in the introduction : 
Is Rayleigh scattering significantly beneficial because of its ability to 
shorten the width of the impulse response? Let us assume that we 
have a slab waveguide with a core-to-cladding index ratio of ni/n 2 
= 1.01. From (59) we obtain in this case 

WaL = 31.2 = 135 dB. (60) 

We may now ask how much loss is associated with a relative decrease 
of the width of the impulse response by a factor 2, or R = 0.5. We see 
from (60) that the amount of scattering loss associated with this 
"improvement" is 

olL = 540 dB. (61) 

This shows that if we are hoping for a reduction in the width of the 
impulse response with the help of Rayleigh scattering, we have to pay 
an intolerably high price in added loss. Since Rayleigh scattering 
losses are known to be quite small, (59) indicates that this mechanism 
does not help to reduce the width of the impulse response under 
ordinary conditions. 

We are thus forced to consider Rayleigh scattering as detrimental 
to light transmission in optical fibers. Fortunately, it is a small effect 
that does not provide prohibitively high losses at visible or infrared 
wavelength. 

It is easy to understand why Rayleigh scattering is not more effec- 
tive in reducing the width of the impulse response. It has been shown 
that a very carefully shaped power spectrum of the function describing 
fiber irregularities is required to reduce the loss penalty for pulse 
width reduction. 2 Rayleigh scattering is particularly poorly suited for 
this purpose since its power spectrum is flat. Only a very small frac- 
tion of the total amount of scattering is used for mode mixing, most 
of it is used for light scattering into radiation modes leading to scatter- 
ing losses. 

Our calculation was based on a slab waveguide model. However, 
the result is expected to be representative of round optical fibers. 
Experience has shown that estimates of the performance of round 
fibers can be obtained from scattering data calculated on the basis of 
a slab waveguide model. 
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